Maxwell-Chern-Simons Hydrodynamics for the Chiral Magnetic Effect 
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The rate of vacuum changing topological solutions of the gluon field, sphalerons, is estimated to 
be large at the typical temperatures of heavy-ion collisions, particularly at the Relativistic Heavy 
Ion Collider. Such windings in the gluon field are expected to produce parity-odd bubbles, which 
cause separation of positively and negatively charged quarks along the axis of the external magnetic 
field. This chiral magnetic effect can be mimicked by Chern-Simons modified electromagnetism. 
Here we present a model of relativistic hydrodynamics including the effects of axial anomalies via 
the Chern-Simons term. 

PACS numbers: 47.75.+I, 11.15.Yc, 12.38.Mh, 31.30.jg 
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Introduction 

The ground state of strong interactions is believed not 
to break CV invariance owing to a theorem by Vafa 
and Witten [l|. However, it has been suggested that 
the Vafa- Witten theorem does not preclude formation 
of metastable V- and CP-odd domains [2|-|5|- These V- 
odd bubbles are expected to occur randomly anywhere in 
space-time in hot quark-gluon plasma Q , and they have 
the property that (E a -B a ) ^ 0, where the fields are color 
electric and magnetic fields 0]. 

P-odd bubbles may yield observable effects on quarks 
in the presence of a strong external (electromagnetic) 
magnetic field because of the winding of the gluon field. 
The necessary setup for observation might be achieved 
in heavy-ion collisions when two nuclei collide with a 
nonzero impact parameter, which results in a strong mag- 
netic field normal to the collision plane. 

In the deconfinement phase of quark-gluon plasma, the 
chiral limit (m quar k = 0) is a reasonable approximation. 
In this limit quarks have definite helicity eigenstates as 
the left-handed and right-handed helicity states do not 
mix in the absence of a mass term in the Lagrangian. 
Accordingly, the momentum and spin will be parallel for 
the right-handed massless quarks and antiparallel for the 
left-handed ones. In the meantime, if the gauge field has 
a nonvanishing winding number, the chirality will not 
be conserved owing to the axial anomaly, and vacuum- 
changing topological solutions will create net chirality. If 
there are more left-handed or right-handed quarks than 
the other, then when they align their spins parallel or 
antiparallel to the external magnetic field according to 
their charge, charge separation will occur in the direc- 
tion of the magnetic field. As a result, a net electric 
current perpendicular to the collision plane will be pro- 
duced. This is called the chiral magnetic effect (CME) 
[8hT^|. Charge separation creates an electric dipole mo- 
ment and it causes CV violation. 

As chirality nonconservation is purely a nonperturba- 
tive effect, any evidence for this effect in the data would 
be a valuable observation of the nontrivial topology of 



Yang-Mills theories. The CME has been verified in lat- 
tice calculations 13- HI]. According to the STAR Collab- 
oration, data from RHIC indicate that charge separation 
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has indeed been observed 



Vacuum of Yang-Mills Theories 

The vacuum of non-Abelian gauge theories are in- 
finitely degenerate; there are infinitely many pure gauge 
configurations of the gauge field. Different vacua are dis- 
tinguished by an integer winding number, defined as 
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where 



e^ uafi G a al3 is the dual field strength. 



No small perturbation around the minimum of a vac- 
uum can change the winding number; perturbation the- 
ory is blind to the many-vacua structure of the non- 
Abelian theory. Instantons are finite solutions of the 
Euclidian action and they are localized in space-time 
with the feature of tunneling between two different vacua. 
Their effect, however, is suppressed exponentially at any 
temperature. Another possibility of the vacuum chang- 
ing solution is the sphaleron. It interpolates from one 
vacuum to another by hopping over the energy barrier, 
whose height is on the order of Aqcd- This is possible at 
high temperatures because there will be enough energy 
to surmount the barrier without tunneling. Although 
sphalerons originally appeared in the electroweak theory, 
their importance also arises in quark-gluon plasma [20j. 



The approximate sphaleron rate is (see Ref. [21| 
references therein) 



r = 



dN 
(Pxdt 
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and 



(2) 



where a s = g 2 /^ and T is the temperature. As 
sphalerons do not tunnel, the rate is not suppressed at 
high temperatures and an appreciable number of transi- 
tions is expected. 
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QED Coupled to QCD with Topological Charge 



Maxwell-Chern-Simons Hydrodynamics 



The charge separation may be examined from a more 
formal viewpoint. The Lagrangian for QED coupled to 
QCD with topological charge is given by [22] 



£qcd+qed — 

$f [*7"(3/. - igAp a - iq f AJ - m f ] 



f 
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g 2 G a ^G^--F^F^, (3) 



where A° is the gluon field, is the electromagnetic 



field, 



and G a ^ u and F^ are the gluon and electromag- 



netic field strength tensors, respectively. Although the 
term F^F^ does not exist in the Lagrangian, it will be 
induced by the topological charge of the gluon field via 
quark loops [23j]. To focus on the electromagnetic sector 
of this theory we can start with the effective Maxwell- 
Chern-Simons Lagrangian 24- 26j|: 



A very strong magnetic field is created in heavy ion 
collisions when two highly charged nuclei collide with a 
nonzero impact parameter. Considering typical RHIC 
parameters, the magnetic field generated normal to the 
collision plane is of the order eB ~ (10 — 100 MeV) 2 dur- 
ing the first moments (r ~ 1 fm/c) of the collision. This 
magnetic field is even stronger than the ones in magne- 
tized neuron stars and magnetars, where typical mag- 
netic fields are about eB « (2 MeV) 2 , or 10 10 T. 

Relativistic hydrodynamics is the prime tool being 
used to describe the time evolution of quark-gluon 
plasma. Let us start with the conservation of energy and 
momentum to derive the hydrodynamic model of quarks 
coupled to Chern-Simons modified electromagnetism. 



d^T*" + en = 0. 



(10) 



Here the energy-momentum tensor of the ideal fluid is 
given by 



£mcs — -^F^F^ 



A^ 



-9F^ V F 



(4) 



where 9 — 9(x,t) will mimic V-odd bubbles, c = 
J2j q 2 e 2 /(2n 2 ) and is the electric current of the 
quarks. The last term in Eq. Q is not a total derivative 
and does not vanish because 9 is a field and is determined 
by fluctuations of the topological charge. The equation 
of motion following from ^ is 



d^ lu = J v - cF av d a 



(5) 



Together with d^F^" = 0, the set of Maxwell-Chern- 
Simons equations can be written as 



(e + P)u"u v - gTP, 



(11) 



where e and P are the local energy density and pressure. 
The fluid four- velocity is u M = (7,7V), and in the rest 
frame of the fluid it reduces to u M = (1, 0, 0, 0). We use 
the natural units h — c = ks = 1 and the metric g^ u — 
diag(+l, — 1, — 1, —1). The energy-momentum tensor of 
the free electromagnetic field is given by 



Qiw = ftp** + Igl^FaftF *. (12) 

Its divergence can be found by using the equation of mo- 
tion ©, 



V x B 



dE 
~dt 



J + c[9B + V9 x E 



(6) 



d„,W = -F vX J x - cF vX F Xa d a t 
Substituting this into Eq. (jTOJ) gives 



(13) 



V • E = p- cV6> • B, 



(7) 



(14) 



dB 
~dt 



0, 



V • B = 0. 



(8) 



(9) 



These reduce to Maxwell's equations when the field 9 is 0. 
The modification of electromagnetism can be attributed 
to the vacuum changing solutions in the full QED x QCD 
theory. 



The electric current of the quarks in the fluid is defined 
as J p = nu 11 , where n is the net electric charge density. 
It is conserved, as can be seen by taking the divergence 
of both sides of §5§ 



= 0. 



(15) 



_ By using the definitions E^ = F^ v u v and B^ 
F^Uv, Eq. (fl4")> can be written as 
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d^T^ = nE v - cE x B x d v L 



(16) 



where E X B\ = —E ■ B. We refer to the right-hand side 
of Eq. (fl6|) as the Lorentz-Chern-Simons force density. 

Now we calculate the entropy production in the system 
by contracting Eq. (fT6"|) with u v . Using the fact that 
u v E v = 0, current conservation, and the normalization 
condition of the four- velocity u^u^ = 1, we find 



e + P\ 1 



n J n 
By using the Gibbs relation, 



- -d^P = -cE x B x u v d v 9. (17) 



dfi±^)=T d | 

n I \nJ n 



\nJ n 



(18) 



Eq. dT7|) can be transformed into 



1 

T 



cE x B x u v d v 9, 



(19) 



where a is the entropy density, and the entropy current is 
= o"ti M . Here we see that the Chern-Simons term can 
increase or decrease the entropy as it can be positive or 
negative. Although it may appear to be a violation of the 
second law of thermodynamics, we point out that entropy 
can decrease locally. However, the change in the entropy 
of the whole system cannot be negative. This situation 
is similar to that of Maxwell's demon, where the demon 
can collect fast-moving and slow-moving gas molecules 
on two separate sides of the container by controlling the 
door in the middle of the container and, consequently, 
decrease the entropy of the gas by driving the system into 
more ordered phases. However, this does not violate the 
second law of thermodynamics. The entropy produced 
by the demon cannot be less than the entropy decrease 
of the gas, which means that the total entropy change of 
the whole system, the gas and the demon, cannot be less 
than 0. In Eq. (fT9|) we isolated the external field 6(x, t) 
and calculated the entropy for the rest of the system. 

Discussion of the entropy also relates to the vacuum 
structure of the SU(N) gauge theories. It is well known 
that the vacuum of non-Abelian gauge theories is in- 
finitely degenerate. Topological solutions, such as instan- 
tons and sphalerons, interpolate between different vacua, 
and what enumerates that is the winding number. As 
positive and negative change in the winding number are 
equally possible, its behavior is completely determined by 
a one-dimensional random walk. Therefore, a net change 
in the winding number can happen in a particular event 
by means of fluctuations. The parallel electric and mag- 
netic fields produced by the term G a ^G a ^ ~ E a ■ B a 
correspond to increase or decrease in the entropy locally; 



however, when infinitely many degenerate vacua are con- 
sidered, the total entropy does not decrease. We note 
that in Ref. [27j an extra term was added to the current 
to cure the negative entropy change problem. Here we 
take a different approach. 

To calculate the entropy of the whole system, we use 
the modified Gibbs relation (see also Ref. |28j): 



e + P 



T<1 ( -) + -dP + -R d9, (20) 
\nJ n n 



where the generalized force Rg is defined as 
Re 



~d9 



7> 



and Q is the thermodynamic potential density. In terms 
of the partition function of the model given in Eq. (QJ , 

Z = J [Vq\[Vq}[DA^]ex-p (j^ dr J d 3 x£j (21) 

where /3 = 1/T, the thermodynamic potential can be 
written as 



n = -p = y\az. 



(22) 



The term -jOF^F^ = c9E ■ B in the Lagrangian, Eq. 

involves fields such as B, which is generated by the 
colliding nuclei, and E, which is induced by the field 
9(x,t). Here we assume that the fields vary so slowly 
in space and time that they can be taken outside the 
integration to give 



Z = cxp(-^c9E x B x 

x J [Dq][Dq][DA'^ V {j\r f d s xC^j (23) 

This assumption is no better or worse than the applica- 
tion of hydrodynamics to heavy-ion collisions, where it 
is necessary to assume that the matter is locally ther- 
malized and that the temperature and flow velocity vary 
slowly in space and time. The functional integration in Z 
involves fluctuations A'^ around the externally imposed 
B field and the induced E field. Hence Rg = —cE x B\. 
Now we rewrite Eq. (|17p as 



mi' 1 ( <),, [ ^-^-) - -d^P + -cE x B x d^9 ) = 0. (24) 



Using Eq. (|20p we deduce that entropy is conserved: 
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= 0. (25) 

In this calculation we used ideal hydrodynamics and ig- 
nored the dissipative terms. It is well-known that the 
inclusion of dissipative terms in T^ u and J p does not 
result in a decrease in entropy. 

In Eq. (fTT)) we calculated the component of the energy- 
momentum conservation equation in the direction of u u . 
Now we find the projection of Eq. (|16|) perpendicular to 
U v . Multiplying both sides of Eq. (fT6|) by the projection 
tensor A" = 8" — u a u Vl we get 

(e + P)u^d^u a - (<S? - u a u u )d u P = nE a 

-cE x B x {&" - u a u v )d v B. (26) 

Defining the comoving derivative as D — u^d 11 , the 
normal derivative as D a = A"d u , and using the defini- 
tion w = e + P, we can write Eqs. (IT71) and (|2^|) as 

De + wd^ = -cE x B x D6 1 (27) 
wDu a - D a P = nE a - cE x B x D a e. (28) 

Equation (|2"8"f is reminiscent of the nonrelativistic Euler's 
equation. 

Conclusion 

We have considered Chern-Simons modified electro- 
magnetism in the context of the Chiral Magnetic Effect. 
We formulated the hydrodynamic model for the purpose 
of incorporating axial anomalies into relativistic hydro- 
dynamics. The equations we found should be solved nu- 
merically by using realistic values for the magnitude of 
the external magnetic produced in heavy ion collisions 
and the size and rate of the sphaleron transitions. This 
work is now in progress. 
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